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The first post-Newtonian (PN) hydrostatic equations for an irrotational fluid are solved for 
an incompressible binary system. The equilibrium configuration of the binary system is given by a 
small deformation from the irrotational Darwin-Riemann ellipsoid which is the solution at Newtonian 
order. It is found that the orbital separation at the innermost stable circular orbit (ISCO) decreases 
when one increases the compactness parameter M*/c 2 a*, in which M» and a* denote the mass and 
the radius of a star, res pectivel y. If we compare the 1PN angular velocity of the binary system at 
the ISCO in units of -J M* / a» with that of Newtonian order, the angular velocity at the ISCO is 
almost the same value as that at Newtonian order when one increases the compactness parameter. 
Also, we do not find the instability point driven by the deformation at 1PN order, where a new 
sequence bifurcates throughout the equilibrium sequence of the binary system until the ISCO. 

We also investigate the validity of an ellipsoidal approximation, in which a 1PN solution is ob- 
tained assuming an ellipsoidal figure and neglecting the deformation. It is found that the ellipsoidal 
approximation gives a fairly accurate result for the total energy, total angular momentum and an- 
gular velocity. However, if we neglect the velocity potential of 1PN order, we tend to overestimate 
the angular velocity at the ISCO regardless of the shape of the star (ellipsoidal figure or deformed 
figure) . 



I. INTRODUCTION 



At the beginning of the next century, we will have powerful instruments for the detection of gravitational waves, 
such as LIGO, § VIRGO, @ GEO @ and TAMA. f§ One of the most promising sources of gravitational waves 
in the sensitive frequency range of these laser interferometers is coalescing binary neutron stars (BNSs). If we have 
accurate theoretical templates of inspiraling phase of BNSs, we can extract various information regarding them from 
gravitational waves, such as their mass and spin, j^] Moreover, when we construct reliable theoretical templates around 
the innermost stable circular orbit (ISCO) of BNSs, we can obtain physical information about the equation of state 
of neutron stars, i.e., the relation between the mass and the radius of a neutron star. || 

Binary neutron stars evolve due to the radiation reaction of gravitational waves, so that they cannot reach equi- 
librium states. However, the timescale of the orbital decay is much longer than the orbital period of BNSs until the 
ISCO. Therefore, we can regard the state of a binary system as a quasi-equilibrium, even if the orbit approaches the 
ISCO. 

From this point of view, several authors have attempted to obtain relativistic quasi-equilibrium configurations of 
BNSs numerically. |m^] However, they assume that the binary system is synchronized. Synchronization (or corotation) 
is not a realistic assumption for the BNSs velocity field just outside the ISCO. This is because the viscosity of a neutron 
star is negligible even near the ISCO, and the velocity field of BNSs becomes irrotational or nearly irrotational . P lGl 

Using the Newtonian theory, Uryu and Eriguchi have recently investigated an irrotational binary system. |ljp 15| | 
However, since BNSs are general relativistic objects, the Newtonian treatment is not sufficient and we need to include 
general relativistic effects. In order to investigate general relativistic effects on the orbital motion of BNSs, equilibrium 
sequences of BNSs composed of Newtonian stars with interaction forces of 2PN order have been studied. |Q On the 
other hand, Wilson, Mathews and Marronetti have computed sequences of irrotational BNSs in the conformally flat 
approximation of general relativity. jlj| They suggest that the central densities of the stars increase when the stars 
decrease their separation and massive neutron stars collapse to black holes prior to merger. This conclusion is in 
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contradicting with that of other works. [fil||8]]l7j] Therefore, it is necessary to solve the quasi-equilibrium problem in 
general relativity for an irrotational binary system. 

In order to obtain accurate theoretical templates of gravitational waves from irrotational BNSs just prior to the 
ISCO, which is our goal, a general relativistic hydrostatic problem with compressible equation of state must be 
numerically solved. Even if we could obtain such numerical solutions, it is necessary to compare them with analytic 
or semi-analytic ones in order to check the validity of the numerical calculation. Lombardi, Rasio and Shapiro have 
semi-analytically studied irrotational BNSs using the energy variational method in the first post-Newtonian (PN) 
approximation of general relativity. However, they do not solve the velocity fields of the binary system at 1PN 
order. They only give the velocity fields of the incompressible fluid at Newtonian order which flows along a plane 
perpendicular to the rotational axis X3. We know that the existence of the X3 component of the velocity field is 
important in the irrotational fluid problem. []lT|-|l3 18 From this fact, we believe the velocity potential of 1PN order 



in the irrotational binary problem is important. 

We have already determined the equilibrium sequence of a corotating binary system. JToj ] However, the formalism 
which we have used in these papers is not applicable for the irrotational case. A formalism for obtaining quasi- 
equilibrium configuration of an irrotational binary system in general relativity has recently been constructed. |2Q-|2^1 
In this formalism, we need to solve only two hydrostatic equations as for the fluid equations. One of them consists 
of the integrated forms of the Euler equation and the other is the Poisson equation for the velocity potential. Thus, 
the formalism seems to be very tractable for computing the equilibrium configuration of realistic irrotational bodies. 
In order to develop the method for solving the irrotational binary problem consistently, the present author and 
collaborators recently investigated an irrotational and incompressible star of 1PN order, |ll| using above mentioned 
formalism. Our method was originally introduced by Chandrasekhar. p5|-|2l| As an extension of this method, we 
study an irrotational and incompressible binary system in the 1PN approximation of general relativity in this paper. 
We assume that each star in the binary system is constructed from an incompressible and homogeneous fluid. This 
assumption presents a great advantage for solving equations analytically. Also, the deformation of the figure of the 
binary system at 1PN order is given by the Lagrangian displacement vectors. 

This paper is organized as follows. In §2, we formulate the method to solve the irrotational binary problem. In 
§3, the Lagrangian displacement vectors for description of the deformation of the binary system are given, and the 
angular velocity of 1PN order is derived from the first tensor virial relation. We give the boundary conditions for 
determining the velocity potential and the deformed figure of 1PN order in §4. The total energy and total angular 
momentum of the binary system are calculated in §5 and numerical results are given in §6. Section 7 is devoted to 
summary and discussion. 

Throughout this paper, c denotes the light velocity and we use units in which G = 1. Latin indices i, j, k, ■ ■ ■ take 
values 1 to 3, and Sij denotes the Kronecker delta. We use and rf y as the quadrupole moment and its trace free 
part, 



±ij - kj - -Sij ^2 !kk (1-1) 

of each star in the binary system. 



fc=i 



II. FORMULATION 



Equilibrium configurations of binary systems with non- uniform velocity fields are obtained by solving the Euler, 
continuity and Poisson equations consistently. Since we consider a binary system composed of incompressible stars, 
all the calculations are carried out analytically, even in the 1PN case. The procedure is as follows. 

(1) We construct a Newtonian equilibrium configuration of the binary system, i.e., the irrotational Darwin- Riemann 
ellipsoid, p9|-|3l| as a non-perturbed state. For simplicity, we consider only the case in which the directions of the 
vorticity vectors and the angular velocity vector lie along a^-axis, and we assume that BNSs are composed of equal 
mass stars. 

(2) The 1PN corrections for the velocity potential and gravitational potentials are obtained from their 1PN Poisson 
equations. 

(3) We calculate the deformation from the Newtonian binary system induced by 1PN gravity using the Lagrangian 
displacement vectors introduced by Chandrasekhar. [29| ] Then, the corrections for the Newtonian quantities due to 
the deformation of the binary system are estimated. 

(4) From the first tensor virial relation at 1PN order, i.e., the force balance equation at 1PN order, we calculate the 
correction for the angular velocity. 19 
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(5) We substitute all the 1PN corrections obtained in (2), (3) and (4) into the 1PN Euler and continuity equations. 
Then, the coefficients of the Lagrangian displacement vectors and the 1PN velocity potential are determined from 
these two equations with the boundary conditions on the stellar surface. 

In this section, we calculate the Newtonian and 1PN terms which we need in the above procedures. We consider the 
equilibrium sequences of BNSs of equal masses (Mi = M 2 — M) whose coordinate separation is R^ For calculational 
simplicity, we use four kinds of coordinate systems which are all rotating frames. The first frame is Xi, in which the 
center of mass of a star (star 1) is located at the origin of the coordinate system, and the other one (star 2) is located 
at (x±, X2, £3) = (-R, 0, 0) (see Fig. 1). The second one is yi, in which the center of mass of star 2 is located at 
the origin, and star 1 is located at (3/1, y 2 , 1/3) = (R, 0, 0). The third one is Xi, in which the origin is located at the 
center of mass of the binary system, and the centers of masses of two stars are located at (Xi, X2, X3) — (-R/2, 0, 0) 
and (-R/2, 0, 0). The fourth one is Y^. This frame is just the same rotating frame as X t , but we label it with a 
different name for convenience. We use X{ for star 1 and Yi for star 2. Then, the relations among the four corotating 
frames are 



/ R \ 
(Xi, X 2 , X 3 ) = [xi + — , x 2 , x 3 j, 



(Yi, Y 2 , Y 3 



yi 



R 



2/2, 2/3 



(2.1) 
(2.2) 



Due to the symmetry, we consider only the equilibrium configuration of star 1 in the following. 



(-R/2, 0, 0) 




Star 2 



Starl (R/2,0,0) 



FIG. 1. Sketch of the Darwin- Riemann ellipsoid. Three coordinate systems are shown. The center of mass of star 1 is located 
at the origin of the corotating frame Xi, and that of star 2 is located at the origin of yi. These two origins of the corotating 
frames are located at (R/2, 0, 0) and (-R/2, 0, 0) in the coordinate system Xi in which the origin is located at the center of 
mass of the binary system. 



The main purpose of this paper is to calculate the 1PN corrections for the angular velocity, the velocity potential, the 
energy, and the angular momentum in the Darwin-Riemann problem. In the calculation, we assume that the parameter 
ao/R, where ao is a typical radius of a star, is small and we use it as the expansion parameter. Accordingly, there 
exist two different types of expansion parameters in this paper. One of them is the post-Newtonian parameter and 
the other is ao/R. 

At Newtonian order, the angular velocity (Odr) becomes |32] (see §111 for derivation) 



"dr. 



2M 



R 5 



0(R- 7 ). 



(2.3) 



^The coordinate condition in this paper is the standard PN one. 
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Thus, in the 1PN approximation, we can expect that the following types of quantities will be the main terms at 1PN 
order: 



MM MM Mai M A Ma l M 

~ W X ' ~ W X W ~ ~W X X ( ' J 

Here we have used the relation -f u ~ Mag. We will derive the four types of terms given above. 

Since we consider an incompressible fluid, the gravitational potentials inside each star are expressed as polynomials 
in the coordinates Xi. For the purpose of obtaining the 1PN corrections given above, we need to take into account 
the coefficients of terms such the gravitational potentials of 1PN order, where < mi, ma, < 5 

and < mi + TO2 + ^3 (= n^t) < 5 up to 0(R~ 5 ) for the case in which m t is odd and up to 0(R~ 3 ) for the case in 
which m t is even. 



A. Hydrostatic and Poisson equations for 1PN irrotational binary systems 

For an irrotational fluid, the relativistic Euler equation can be integrated, and in the 1PN case, it is written as (22 
constant =--U+\ ^(d fc </> N ) 2 - ^ fe 



—u + \u 2 + X - 



k 

1 (P 



2\p 



+ ^(^n)(^ pn ) - J2e k d k 4> PN -J2PA<j) N 



(2.5) 



where we have assumed that the fluid is incompressible, i.e., p = const. In Eq. (|2.5|), P, p, £ k , </>n, <?W, U, X, and 
Pk denote the pressure, the density, the velocity field of the orbital motion (spatial component of the timelike Killing 
vector), the Newtonian and 1PN velocity potentials, and the last three terms are the Newtonian and 1PN potentials, 
which are derived by solving the Poisson equations as 



AU = -Airp 
AX = Airp 



2{/+-+2V(^ N ) : 



AP k = -inpd k (t) N , 

A x = ^ P Y^{dk^)X k . 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



Here (3k is expressed as 



$k = - 



\Pk + \(d kX + Y, xld kPi 



(2.10) 



We note that using the gravitational potentials, the spacetime line element to 1PN order can be written as 

ds 2 = -a 2 c 2 dt 2 + -| hdx l dt + (l + ^2") dx%dx ^ 



(2.11) 



where 



U 1 fU 2 \ 



-6\ 



(2.12) 



A characteristic feature of the irrotational fluid is that the continuity equation reduces to a Poisson-type equation 
for a velocity potential <f>, 0,0] and in the 1PN incompressible case, it becomes 



(2.13) 
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where 



Ci = -t + d i( fa + 



l r ™/i 



(2 E( a ^ N ) 2 + 3C/ ) - ~ & + 9 ^ PN 

k ? 



(2.14) 



In the following subsections, we obtain the terms appearing in Eq. (2.5) by solving equations separately. 



B. Newtonian terms 



Each gravitational potential is composed of two parts. One of them is the contribution from star 1 and the other is 
from star 2. In the following, we denote the former part as ^P 1 ^ 1 and the latter one as \t 2 , where *f? denotes one of 
the potentials. We also define v? 1 ^ 2 and vf 2 ~* 2 as the contribution from star 1 to 2 and star 2 to itself, respectively. 

In analogy to previous studies, ]29|-|3l[ the configuration of each star in the binary system at Newtonian order is 
assumed to be an ellipsoidal figure of its axial length eti, 02, and 03. In this case, the solution of the Poisson equation 
for the Newtonian gravitational potential 



AU 



-47rp, 



(2.15) 



is written as U = U 1 



U 



U 2 



1— >i 



L , where 



tt P {A -J2 A 



M 
R 



1 

'-I 



0t. 



O ,-•-> 2 2 ,-vt 2 



-2a:? + Zx x {x\ + £§) 



2i? 2 2R? 
8a;f + 3x| + Zx% - 2\x\(x\ + x\) + e^! 1 



8i? 4 

3^ u / 3jbi 12a; 2 -5a; 2 - 5a; 2 \ 3 . a _ 2 . 

+ ^pg (8/1111 + 3/2222 + 3/3333 — 24/n22 — 24/n33 + 6/2233) 

+o(ir 6 ), 



(2.16) 



(2.17) 



and 



d xpx i x 



M 

35 S ' ' 



a 2 a 2 (l + 2%). 



(2.18) 



Here, Aij... are index symbols introduced by Chandrasekhar |2{| and Aq = J2i Aiaf is calculated from |2£ 

A Q = aia 2 a 3 



y (of + u)(a| + u)(a| + u) 
a/(1 + t)(a| + £)(a 2 + *)' 



(2.19) 



where 02 = 02/01 and 03 = 03/01. Note that U 2 ^ 2 and t/ 1 ^ 2 are obtained by changing xi to — (x\ + R) in J7 1 
and U 2 ^ 1 , respectively. 

The pressure at Newtonian order is written as 



p = Po(i-£ 



9 1 1 

1 a r 



(2.20) 



where Pq denotes the pressure at the center at the star and is calculated from the scalar virial relation as 

- = \*pM - ^11 - \F a {al - a 2 2 )Q 2 DR + 0(iT 5 ). 



(2.21) 
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The axial ratios 012 and a 3 are determined from plfl 



D TUT „J. 

a\Ai = l± + ^ + -^ R F a (2 - F a ) + 0(R-% 

^^-^ + 0(rS) ' (2 - 22) 



We can determine the equilibrium sequence of the irrotational Darwin-Riemann ellipsoid from Eqs. ( [2.3] ) and ( 2.22 ). 
The velocity field of Newtonian order of star 1 in the inertial frame is written as Vi = 9i</>N, where 

2 r 

Vl =£ 1 +ui = -( a 2 llR 2 + l)nx 2 , 

V2=£ 2 +U 2 = ( ° 2 ^ ? + lWl + ^fi 

V a\ + a 2 / 2 
v 3 = 0. (2.23) 

The quantity m is the internal velocity of the star in the corotating frame xt, and f2 denotes the angular velocity of 
the binary system. £i and Ui are given by 

ii= (-Clx 2 , n(aJi + |), 0), (2.24) 

Ui = (—Ax2, ~ —Axi, 
Vq.2 a\ 

= (-4Kvx 2 , 4h^i,0), (2.25) 

where A is the angular velocity of the internal motion. Jr is defined as C/^i where £ = (rotu)3 denotes the vorticity 
in the corotating frame. For the irrotational case, /r becomes —2. Thus, the velocity potential at Newtonian order 
becomes 

R 

0n = F a ilx 1 x 2 + —Slx2 (2.26) 



and the velocity field is written as 



where we define 



(F a n X2 , F a n Xl + ^n, 0), (2.27) 



F - s 3+3- < 2 - 28 ' 



C. 1PN terms 

1. X 

For computational convenience, we separate X into two parts as 

X = X + X V , (2.29) 
where X and X v are derived from the Poisson-like equations 

/ 3P\ 

AX = 4ttp[2U+ — J, (2.30) 
AX V =8TrpJ2(dkM 2 - (2.31) 



G 



As in the case of U, the 1PN potentials are divided into two parts as Xq = X^ 1 + X 2 ^ 1 , and we consider them 
separately. 

The contribution from star 1 is derived from the Poisson-like equation 



Axr 1 



3Po 



and the solution is 



A 



-aoll 1 ^ 1 + a 1 D 1 + ^\/ ; L> 



M 

„ - ^(2D n - D 22 - D 33 ) 



where 



-^(2Diii - 3D 122 - 3D 133 ) + 0(R- 5 ), 



„ „ 3P 2M 3^ii c , 

a = 2irpA + — + — + — £ + 0(R~ 5 ), 



Oil 



P 

2M 9^ii 
9 A x 3P ° 



i? 



0(iT 



I? 3 



Here, are the solutions of equations 



AD, 



-AnpxiXj 



and the solutions are written as |29j 

i 

1 (An — S ' Anixf^x 



Da = irp 



i 



l l m 

D i:j = npafa 2 (A tj - 2J Aijix 2 )xiXj, (i ^ j) 



2 2 s ! 
ilmXi X m j 



•f — Aiiiixf)xf 

I 

3 ^ ^ ^ 



Dj 



Ttp 



ajctj(Aijj — ^ Aijjixf)xiX 2 



l m 
2 



/ m 



where Bijf.,,, are index symbols defined by Chandrasekhar. |29| 

The contribution from star 2 for is calculated from the equation 



AX^ 1 = 4np 



(2,pA + ±f)-j:{2npA l + ^)y? + 2U 



l->2 



where j/i = — (ajj + -R), 2/2 = %2, and 2/3 = x 3 . The solution is written 



Y 2—>1 



1 

— {201^ - AY 1 - Dtt 1 ) + 0(R-% 



(2.32) 



(2.33) 

(2.34) 
(2.35) 
(2.36) 

(2.37) 
(2.38) 

(2.39) 
(2.40) 

(2.41) 

(2.42) 



(2.43) 



(2.44) 
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where D 2 ^ 1 are calculated from the same equations as in the case of A,..., i.e. 



The solutions are 



where 



D 



D 



° l ~ R? 



ADf~! = \ -/,!,,,,, • • • . 



1 ~ IT + 2P? ~ X 2 - A) - ^3(^1 _ ix ^ X 2 - ?> X l X l) 



+ 2^4(2/1111 - 3/ii 22 - 3/1133) (l ~^f) + 0{R~ 



R 



ti {l-^)+0(R-% (i^l) 



2/? 2 



2R 3 



+ ^(i-f) + o(«- s ). 



^.„f(i--2L) + o(« 



-6\ 



(Ml) 



2->l _ -'Hjj 
1« fi 2 



l-^+2^(2x?-x»-xI)] +0(i?- 5 ), 



2^1 _ ^22ji 



^2jj - 



/? 2 

2->l _ -^33jj 
3M ^2 



■X2 3xix 2 " 
x 3 3xix 3 



+ o(/r 5 ), 

+ 0(R~ 5 ), 



iill — -I-iill 



o y 1 1an • 



(2.45) 



(2.46) 
(2.47) 

(2.48) 
(2.49) 
(2.50) 
(2.51) 
(2.52) 

(2.53) 



As in the case of Xq, it is convenient to separate X v as X^ 1 + X 2 ~^ 1 . Here, X^ 1 is derived from the Poisson-like 
equation 



AX^ 1 = 8ir P n 2 F 2 {x{ + x 2 2 ) + F a Rx\ 



R 2 



and the solution is written 



The equation for X 2 ^ 1 is 



-i->i 



-2^ 



F 2 (D n + D 22 ) + F a RD 1 + ^-U 1 ^ 1 



AX^ = Snpn 2 F 2 (yf + y 2 2 ) - F a Ry 1 



R 2 



Then, the solution is found to be 

x 2 r x = -2n 2 



F^Djr 1 + Dlr 1 ) - FaRD 2 - 1 + ^-U 2 ^ ! 



R 2 
4 



(2.54) 



(2.55) 



(2.56) 



(2.57) 



2. p k 

Substituting Vi of Newtonian order into the equations for Pj and x, we immediately find the solutions 
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p^ 1 = F a nD 2 , 



P i-i 

r 2 

pl^l 
p2^1 



= 0, 



F a nD 2 2 -\ 



p2-,l _ 
r 2 — 



x 1 - 91 = -n 



x 2 - 1 = -n 



2F a D 12 + ^(F a + l)D 2 



2F a D 2 ~ 



Using the solutions of Pk and % we have derived above, we obtain 



where 



Pi" 1 = -^np£lx2 



a\A x - 7a\A 2 - 2a\a\A 12 

(5a^A 12 — ia\An + Qa 1 a 2 Au 2 )x 1 



+{7a%A 22 - a\A 12 + 2a\a\A X22 )x\ 
+ (7a 2 2 A 23 - a 2 A 13 + 2a 2 a 2 A 123 ); 



R 



— vr pQ(A 1 - ajA 12 ) Xl x 2 , 
F r 

Pl^ 1 = y?rpOa;i a\A 2 - 7a\A\ - 2a\a\A 12 

+{7a\A 11 - a\A 12 + 2a 2 1 ajA 112 )x 2 1 
+ {ba\A 12 - ia\A 22 + §a\a\A 122 )x\ 

+ {7a 2 1 A 13 - a\A 23 + 2a\a 2 2 A X23 )x 2 3 

R r 

+ —npfl -7A - a\A 2 + (7A X + a\A 12 )x\ + (5A 2 + Za 2 2 A 22 )x 2 2 
+ (7A a + a 2 2 A 23 )x 2 3 , 
Pi"* 1 = -^Tr P nx 2 x 3 2F a (a 2 1 A 13 + a\A 23 - 2a 2 1 a 2 2 A 123 )x 1 
+R{A 3 - a 2 2 A 23 )^ 

-Mx 2 



0i^ = 



- !§ (2F„7ii + WF a I 22 + -loo 



\ M 
) ~ R2 



^22 



pr 1 = 4 



2R 

+^ X2 (2xl-xl-xt) + 0(R- 4 ) 

+x i{-™ + U UFJ "-H 2F " + l 

7M n „ „, M 



X\X 2 



637- 



)} 

, 2 2 2 M 2 7M 2 o 2 o 2^ 

-f|^ + o(ir 4 ) 



~ 2 ^! _ MO 

^ 3 "li? 2 " 



2i? 3 



i-^ + o(i?- 2 ) 



(2.58) 

(2.59) 

(2.60) 
(2.61) 

(2.62) 

(2.63) 

(2.64) 

(2.65) 
(2.66) 



(2.67) 



(2.68) 



(2.69) 



(2.70) 



(2.71) 
(2.72) 
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Then, the divergence of (3k is 



where 



^2 dkPl^ 1 = Znpfl 2F a {a\ + a\)A X2 x x x 2 + RA 2 x 2 

k 
k 

Here we have used some relations among index symbols. [E9| 



(2.73) 

(2.74) 
(2.75) 



3. 0PN 

An explicit form of the Poisson-like equation for <^>pn can be written as 

A0pn = 9, 



''*( l^ d M 2 + 3C/) + -difa + ft 



2P a? - 



l l u 2 



a|(a?-a|) n2 | P 



pal 



■1-2 



+n x o(i?- 4 ). 

The solution of this equation can be written up to biquadratic terms in xi a^] 

0pn = (p + qx\ + rx\ + sx\)x\x 2 + (e + fx\ + gx\ + hx\)x 2 + const, 
where q, r, s, /, g, and h satisfy the conditions 



3g + 3r + s = — 



f + 3g + h = -^n 



Pq of - a| 

2 2 



0, 



a 2( a l - a l) 



{al + alY 



n 2 



(2.76) 
(2.77) 

(2.78) 
(2.79) 



On the right-hand side of Eq. ( 2.77 ), coefficients p, q, r, and s are concerned with the star itself. We can find these 
coefficients in a previous paper [18[ in which we investigated irrotational and incompressible stars. On the other hand, 
the coefficients e, /, g, and h are associated with the binary motion. 



D. Collection 



Substituting the terms derived in the previous subsections into Eq. fl2.5| ), we obtain 
P 



g = — u-su 
p 



n 2 

"DR 



T sn- 



r R i 

\F a {2 - F a )x\ - F a (2 + F a )x\ + Rxt + —J 



1 r 



c 

= const, 



to + E + E ^^^^m + ( K ° + X! k ^ + E Kimx2 i x ™) 



Ki 



Kr 



(2.80) 



"'"Although we can add higher order terms which satisfy A<j!>pn = 0, we neglect them for simplicity because we consider only 
biquadratic deformation of ellipsoids in this paper. 
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where S^dr denotes the angular velocity of an irrotational Darwin-Ricmann ellipsoid, and <50 denotes the 1PN cor- 
rection for the angular velocity. 7^ and Ky are expressed as 



7o = 7TP 



P 



«o , Po 



Tip TTP* 



M 3^1 



R 2R 3 



aoA 



2i? 3 



M 
4R 3 



^(2a;Bi - a 2 2 B 2 ~ a 2 P 3 ) - -fi 2 i£(a?Bi + a 2 P 2 ) 



1 



P 
i 

3^ii 



tfR 2 (^-alA 2 )+0(R- 5 ) 



M 

+ 2IlA M R 2 

n 2 



hi R n4 



UMR+ i{l6F 2 (7n + / 22 ) - 30F a /n + (2F a + l)/ 22 + 21^ n }] 



+0(R- 6 ) 

71 = 



/a P ,\^7 2/2. 1 ,-, \ 1 2 ,-> 

~ W + ^ " ° ) & + m<l1 11 2 Sn J " S^* 12 



-^3 (4afAn - 2a 2 Bn + a 2 Pi 2 + a 2 Pi 3 ) 



-2°{-( a 



2^12 _ 



irp 2 a\ 



+ F, 



+Aa\A 11 - 2a 2 a 2 A 12 - 2a 2 5 n - 2a|P 12 + — % 

Ttp z 



)} 



+(F a -m P + 0(R- 4 ), 



72 = 



(Q0 _ » Mo)jl2 + ^ (_ 4o + 4A 2) - ^(a\A, - A)(| + gl) 

+ ( 



7Tp 7T/3^ 



M 
2~R3 



\^] X a\B X2 + T] 2 a 2 2 {a\A 22 - 7^22) 



M 



-lo 2 {^(3a 2 A 22 - 2A 2 - -^1 + F H 3A ° + a 2 A, - 7a 2 A 2 
Z 1 4 \ ix p a 9 / v 

)} 



7r,cra 2 < 

+4a|A 22 - 2a?a|Ai 2 - 2a 2 5i 2 - 2a|P 22 
+ (F a + l)f!p + 0(ir 4 ), 
73 = np 



7T/9 2 



„ x „ Po , „ 1 / 2 „ Pa\(M 3^ii \ 

(ao - irMo)* + -^(Ao + a 2 A 3 ) - ^ (a 2 A 3 - — 2 ) (- + — ) 



+ 



^0 



M 
2~R3 



-m a i B i3 



1 



-77 2 a 2 B 



21! 



+r/3a 2 (a 2 A 3 3 - ^33) -^ElA 



(2.81) 



(2.82) 



(2.83) 
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723 = 7T/0 





<M_ _ 




[-7 




a 2 J 



M 



T^p 2 a\ irp 2 a 2 . 



z \ z J 

-7730^0^233 - 2 fi 233) 

~2R? ^°2^ 223 2 fl 3^233 + 2a^Bi23 — «2-^223 — a 3^233) 

-ifi 2 F a 2 (-3 A3 + 7a 2 A 23 - a 2 As + 2a 2 a 2 A 123 - 4a|A 223 

+2a 2 P 123 + 2a 2 P 223 - —5-^) + 0(iT 
+ (F a + l)Qs + 0(R- 6 ) 1 

A , r> 1 

733 = 7T/0 

1 



-Trp— I — - 



k = irp 



/ 7 

1*a\- M ( P ° 

I ' ^ \np- 2 fl3 3 J 2a 2 i?3 Up 2 " V 
+ i?7ia 2 Bi 33 + i?7 2 a 2 B 23 3 - 773a 2 (a 2 A 333 - ^£333) 

— ^3 (4a 3 A 333 + 2a 2 _Bi33 — a 2 _B 233 — aaPaaa) 

-io 2 P 2 (a 2 P 133 + a 2 B 23 3) + 0(R 
/ P a a ,\(M 9^ii 



M 9^ii \ 2 yr^ni (hi . 9^hii\ 

+|^4 a ? (2a 2 Bn - a 2 Si2 - a 2 3 B 13 ) 
1. 



-io 2 P a p(-A + a?A - 2a\a 2 2 A 12 + ^) + 0(p- 6 ) 
4 V " np A / 

M 2 3M 2ai 

2(13P a - 7)M - -^{4P a (llP a - 15)/ n + (12P 2 + 2P a + 3)J : 
+3(21-13F a )^ n } 



22 



_ l fi 4 P 3 P a + (P a - l)fie + 0{R- 7 ), 
Id 



/ P . a . \ M 1 / P, 
l^ 2 " + ^" A °^"-(- 



\ M 1 / P 2 \ / M 9^ii 



k 2 = np 



+ ^a 2 1 (4aiA 111 - 6a 2 P m + 3ajB 112 + 3ajB 113 ) 
--0 2 PP a (A - fl?4 u + 2a 2 a 2 A 112 - — ^ J + 0(iT' 

~W + 4^° 2(6F ' - 13Fa + ?) - \ Q4RF « + {Fa 1)0/ 

+o(p- 7 ), 



P a , \3M 1 / P 2 \/M W 
V^ 2 + ^ - 2P? " 4 \^ ~ a2A V VP 2 ~2~fi 



(2.89) 



(2.90) 



(2.91) 



(2.92) 
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— 2^4 a i(2 fl 2 y li22 + 2a^i3ii2 — a\Bi22 — a 3 Bi 23 ) 

1 / IP 

+ -n 2 RF a (2A 1 +A 2 + 2,a\A 12 - 2a 2 2 A l2 - Qa\a 2 2 A 122 + —± 

4 V 7Tj0^t 

+0(iT 6 ) 
9 M 2 1 M 1 

~w + H^ 2(4F * + 5Fa - 6) - l nRF « + 2{Fa + 1)0/ 

h3(F a - l)Slg + 0{R- 7 



K 3 = Tip 



3 ^ 



V irp 2 irp 
li3 

1 



tt« °/2i? 4 a|V7rp 2 3 7VR 2 2R 4 



27rpi? 4 , 

3Af 

^^-a 1 (2a 3 Ai33 + 2a 1 S 113 - a 2 B 123 — a 3 B 133 ) 



n 2 RF a (A 3 - a 2 A 13 + 2a\a\A 123 - 



np a 3 1 



9M 2 M 

--p- + ^ 2 ( 13 ^ - 21 ) + - i)nA + o(r- 7 ), 



Mirp 



i? 4 



p o , 1 

5— 2 + A i 

np z af 



K12 = 



M7rp 

«13 = 



a 2 

-^(8a^mi - 6a?5 mi + 3a|S 1112 + 3a^ m; 
AfTrpr P /2 3x 

-a 2 (4a 4 ^m2 - §a\A 1122 - 6afB 1112 + Za 2 2 B 1122 + 303^1123) 



2i? 4 





/ 2 




L TTp 2 




-A) 



2i? 4 L rrp 2 \a 2 a 2 J 

-a 2 (4a 4 ^iii 3 - 603^1133 - 6afBn 13 + 3ajBn 23 + 303^1133) 
_ 3Mnpr P 

K22 -^r^[^4' A2 



a 2 

~-^{^ a 2 A\ 222 + 2a 1 Bn22 ™ a 2 B\ 222 — a 3 Bi223) 



3Mtt P 









Lttp 2 


^2 " 


a 2 / 



«33 - 



+a\(2a\A 1223 + 2a 4 A 1233 + 2a 2 J B 1123 - 02^1223 - «3 B 
Po 



3^1233, 



7rp 2 a 2 



A3 



(2.93) 



a 2 

+ "^( 4a 3^i333 + 2a 1 Bn33 - a 2 Bi 233 - a 3 B 1333 ) 

The quantity SU denotes the gravitational potential induced by the deformation of the binary syster 
form of SU is given in Appendix [a| . 



(2.94) 

(2.95) 
(2.96) 
(2.97) 

(2.98) 
(2.99) 

(2.100) 
The explicit 



III. THE POST-NEWTONIAN ANGULAR VELOCITY 



In this section, the 1PN correction for the orbital angular velocity at 1PN order is derived using the first tensor 
virial (TV) equation. The first TV relation is derived from 
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dP 
d 3 x— = 0. 
ox i 



(3.1) 



Substituting Eq. ( (2.80) ) into Eq. (3.1), we have Eq. (2.3) for Newtonian order. At 1PN order, the explicit form 
becomes 







MR 



5n 2 +6 J d 3 xp^- + Q 2 DR F a (2-F a ) J d 3 xp£ 1 



— [KqM + 3ki/h + K2/22 + «3^33 + 5/tllJmi + 3Ki2-/il22 + 3fti3-Zii33 
+ ^22^2222 + ^23/2233 + K 33^3333)) 

where £1 is the x\ component of the Lagrangian displacement vectors £i. 



(3.2) 



A. Deformation of the figure 



The density profile of an incompressible and homogeneous sphere at 1PN order is the same as that at Newto- 
nian order. However, the density profile of an ellipsoid at 1PN order is different from that at Newtonian order. 
Chandrasekhar's method to obtain the 1PN correction to the Newtonian figure is as follows. |27j] First, the ellipsoidal 
figures at Newtonian order are constructed. Next, the 1PN effect is regarded as a small perturbation to the Newtonian 
configuration, and the deformation from the Newtonian ellipsoid is calculated by using the Lagrangian displacement 
vectors. Finally, by solving equations for the Lagrangian displacement and calculating the correction to U by the 
deformation, the equilibrium configuration at 1PN order is obtained. In this paper, we follow this method. 



In choosing the Lagrangian displacement vectors £^ , we require them to be divergent free J\ ^fc£fe 
incompressibility. To obtain Sil 2 up to 0(R~ 3 ) x fi^, we write the Lagrangian displacement vectors as 



to) 



to) 

k > 



0, due to 



(3.3) 



where 



,,(31) 
.(32) 
.(33) 
t(0) 
.(21) 
.(22) 



{xi, 0, -X3), 
(0, x 2 , ~x 3 ), 

(^1, - XlX2: o), 



x\x 3 



X 3 X\ , 



3 3 



(i, o, 
(5*?, 0, 

(0, xix 2 , 



XiX 3 ^j, 
XlX 3 ). 



(3.4) 
(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 
(3.11) 



Here, we consider only up to the biquadratic deformation so that the Lagran gian dis pla cement vectors of higher 
order functions in Xi can be neglected. The Lagrangian displacement vectors (^4) - (|3.8| ) are concerned with the 
defo rma tion of the star, i.e., coefficients of the velocity potential of the 1PN order p, q, r, and s. On the other hand, 
(3.9) - ( 3.11 ) are associated with e, /, g, and h. 

In order to obtain the orb ital angular velocity of 1PN order, we must calculate the second and third terms on the 
right-hand side of Eq. (3.2), which come from the displacement of the fluid element by the 1PN correction. The 
second term can be evaluated as follows: First, the contribution from star 1 is zero because 



S / d 3 xp 



dU 1 



dx\ 



d 3 xp(x) J d 3 x'p{x l ) 
= - I d 3 xp(x)5 [ d 3 x'p{-x!) 



Xi 



I*- 

Xi 



,/|3 



,/|3 
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d 3 x'p(x')6 f rf 3 xp(x) ^- ^ 3 
J |x x I 



0. 



The contribution from star 2 can be separated into two parts as 

T 2^i r _A d 2 U 2 ^ 1 



5 I d A xp 



8U 2 



dxi 



i=l 



dxidxi 



d xp 



dSU 



2-»l 



dxi 



(3.12) 



(3.13) 



The first term on the right-hand side of Eq. (3.13) denotes the force to which the displaced element of star 1 is 
subject from the nondisplaced potential of star 2. On the other hand, the second term denotes the force to which 
the nondisplaced element of s tar 1 is subject from the displaced potential of star 2. For the Lagrangian displacement 
vectors given in Eqs. (3.4) - (3.11), the first and the second terms on the right-hand side of Eq. ( 3.1 3| ) are equal, and 
we can calculate as 



au 



S J d 3 xp— = 2 J d 3 xp^2c 



2rr2->l 



d 2 U 



dxidxi 



i=l 

d A xp— 

OXi 

M r 

--^{3Sii(2/ u + l33 

+S32 (3/2233 — J2222 
„ ,,/M 18^11 \ 

M t_ 9^1111 , 9^ii/ii , 12/1133 



35*12 (^33 — -^22 
533(13333 



I" 53i(2Jiin 
6/ll33)} 



3/ii22) 



R 2 



MR 2 



R 2 



9 12M (I T \ 

-J22 „r 1^1133 — Jll22j 

R° 



Also, the third term on the right-hand side of Eq. (3.2) is written as 

n 2 UR F a (2 - F a ) J d 3 xp£i = ~I% K F a (2 - F a )(MS Q + 5 21 / n ). 



(3.14) 



(3.15) 



B. The post-Newtonian orbital angular velocity 



The angular velocity at 1PN order is calculated up to 0{R ) from Eq. (3.2) as 



/? 5 



R 
25 2 i 



35n(2/n + J33) + 35i 2 (/ 33 - h-z) + 5 3 i(2/ mi + 3J U22 ) 

+ 532(3/2233 — /2222) — 533(/3333 + 6/1133) 
^DR^a(2 - F a ) 



2S \2M , 36^11 , 1 . 



R 3 



R 5 



2/n IS/ 1 !! 2 IS/ 1 !!!! 24/n 33 1 - 2 
01 + ~ + ^ 5 + ^^DR^a(2 - F Q )ai 



/? 3 5/? 5 



R 

48 

"~Dfi 522(/ll33 — /1122) 



R 5 



R e 
+2np 



f a o 4 . 2 „ \/M 9^ii \ 2ai 2 , 



9^11 / Po , , 
^0 



2/? 5 V7T/9 2 7T/9 



9/^11 

TrpV/? 3 ' 2/? 5 



1G 



-2^h^ E ^ + ^dr^(2«^ 2 + alA 3 ) 

+2(2afA 1 -4A 2 -a^ 3 )} 

2M 2 2M 4ai 
— ~ "R6"*- 8711 ~ 3/22 ~ 3/33 ' + IF 111 

-l04 Ri? 2 Fa | 1 + _4_ (3a 2 + a 2 )F 2j 
-^DR.{(13^a-7)M 

--^ (2(20F a 2 - 41F a + 21)/ n + 3(4F 2 + 5F a - 6)/ 22 

+(13F -21)/3s)} 
, 2£!dr / 3 9 3 9 1 9 \ 4a, , 

+ - !)-^( e + + 5 4? + 5«3^j + gjf + 1)<W 

+0(iT 7 ), 

3Sii(2/n + J33) + 3Si 2 (/ 33 - J 22 ) + S31 (2/iin + 3/n 22 ) 



4 

#5 



+ S32(3/2233 — /2222) — <S , 33(/33 3 3 + 6/1133) 

2S 



R 

2S 2 i 



i? 3 



i? 5 



i? 

48 



2/n 18/n 2 
+ _„ c Qi 



1 



i? 3 



5i? 5 



+ + w nuM2 Fa)a \ 



«S' 22 (/ll33 — /1122) 



N 20 DR / 3, 3, 1 9 \ 4a 2 , , 
+ (F a - + + 5«2.9 + 5«3^J + j£(F a + l)fi DR / 

+ ^(4 + F a )A + f^(2 - 3F a ) + ^(11 + 3i^ / n 



M 



2(-3F 3 + 12F 2 - 77F a + 54) Jn + (-10F 3 + 63F a - 43)/ 2 
+ (55F a - 46)/ 33 l + 0(R~ 7 ). 



(3.16) 



(3.17) 



IV. BOUNDARY CONDITIONS 

In this section, we derive equations to determine the coefficients of the 1PN velocity potential and the deformation 
of the binary system, p, q, r, s, e, /, g, h, and Sij, from the boundary conditions on the stellar surface. The 
stellar surfaces of the Darwin-Ricmann ellipsoid and its deformed figure are expressed as Sur{x) — and S(x) = 0, 
respectively, where 



Sdr(z)=£3t-1, 
i a i 

x 2 



(4.1) 
(4.2) 



Sbr(x) , 5n I -o 2 + S\2 I —o 9 + ^31 I ^—9 2~ 

c 2 1 \aj (13/ Va 2 a 3 / \3a( a 2 / 
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4 5 v + 



2 2 



Sr) +533 ( 



2 2 



3a§ 



S21 



+ 



S 22 ( 



2 



2a 2 
a§ ) 



4 1 



The boundary conditions for the continuity equation and integrated Euler equation are, respectively, 

on 5 = 0, 
on 5 = 0. 



(B) - = 
P 



(4.3) 

(4.4) 
(4.5) 



The condition (A) comes from the constraint that the no rma l com pone nt of the velocity on the surface vanishes. The 
condition (B) determines the stellar surface. Equations (|f.4|) and (4.5) must hold to order c~ 2 . 



A. Condition (A) 



Substituting Ci given by Eq. ( 2.14 ) into Eq. (4^), we obtain the equation at 1PN order 

= - — (a Q + aix x + a 2 x + a 3 x 3 ) 



+ 



2VLx 2 



(A + \ix\ + \ 2 x\ + A3X3), 



(4.6) 



where Uk {k = 0,1,2,3) are functions of a^, p, q, r, s, and Sy, and also Xk (k = 0,1,2,3) are functions of a^, e, /, 
g, h, and $y. Equation (4.6) must be satisfied on Sdr- Then, we obtain the following six equations (three equations 
for p, g, r, s, and Sij, and three equations for e, /, <?, h, and Sij): 
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We can derive the following two equations which do not contain Sn, S\ 2 and p: 



0. 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



,(A) 2 



a z 2 = 0, 



E[ A) a 2 - Z^a 2 = 0. 



Thus, Eqs. (4.13) and (4.14) depend only on the coefficients q, r, s, S31, S 32 and S33, i.e., the 
biquadratic deformation. Also, we can construct two equations which do not contain e and Sq as 
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(4.14) 

coefficients for 
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0, 
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(4.16) 
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Thus, Eqs. (4-15) and (4.16) depend only on the coefficients /, g, h, S21 and S22, i.e., the coefficients for cubic 
deformation. The latter feature is important in determining coefficients. 



B. Condition (B) 

Substituting SU, which is given in Appendix |A|, and the surface equation S = into Eq. (2.8C), we obtain 
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Since (P/p)s must vanish on the stellar surface, we have eight conditions, 
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We can see that the equations for determining the coefficients of the 1PN velocity potential and the deformati on o f 
the b inary sys tem are separated into two combinations. One of them is constructed from Eqs. ( 2.78j ), ( p~7| ) - (4.E) 
and (pL3l| ) - flog) , and deter mines the coefficients p, q, r, s, Su, S12, S31, 6*32, and S33. These coefficients are 
associated with the star itself, and the tripla ne sym met ric pa rt of t he L agra ngian displacement vectors. The other 
is the combination composed of Eqs. ( [2.79D , ( 4.10 ) - ( 4.12| ) and ( 4.36j ) - ( 4.38 ). These equations determine the 
coefficients e, /, g, h, So, S21, and 822- These coefficients are associated with the binary motion, and the 7r-rotation 
symmetric part of the Lagrangian disp lacement vectors around the X3-axis. 

In the case o f an isolat ed sta r, [|18| we have two identical equations for Eqs. ( 4.31 ) - ( 4.35) whe n we su bstitute 
Sfl 2 into Eqs. (4.34) and ( 4.35| ). However, for the case of a bin ary s ystem, the five equations (4.31) - (4.35) are all 
independent up to 0{R~ b ) even if we substitute 5£l 2 into Eqs. ( 4.34 ) and ( 4.35| ). 

Contrastingly, when w e sub stit ute Sf l 2 into Eq. (O6), it is found that So vanishes. This implies that only six of 
s even conditions ( 2. 79] ), ( 4.10 ) - ( 4.12 ) and ( 4.36 ) - ( 4.38 ) are available. Here, we choose equations excluding Eq. 
( p6| ). 

There remains the problem of determining the remaining two coefficients e and So because we have only one equation 
(4.10) (or (4.11) or (|4 . 1 2|) ) for these two variables. There is one degree of freedom. This fact follows from point made 
by Chandrasekhar |26|J27| ] and Bardecn |33| for an isolated star that there is no unique definition of the 1PN solution 
as the counterpart of a Newtonian binary solution. 

Since the condition for connecting a 1PN solution to a Newtonian solution can be arbitrarily chosen, in this paper, 
we simply give the condition as 



S = 0, 



(4.39) 



i.e., we fix the 1PN correction to the orbital separation of the binary system. Then, we have two equations, ( 4.10| ) 
and (4.39), for solving for e and So- 
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Note that when we calculate Eqs. (4.34) and (4.35), we neglect terms such as 
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F a x (e, /, g, h) 



(4.40) 



in order to separate variables into two groups. We take terms up to 0(R~ 3 ) in Eqs. (4.34) and ( 4.35| ). In this 
calculation, we only count the orde r of l/R and ignore the order of the implicit terms, such as F a which is 0(R~ 3 ). 
However, we regard the term (4.4C) as 0(R~ 4 ). Even if we include this term, the physical values, i.e., the e nergy, the 
angular momentum and the angular velocity, do not change. § Note that when we include the term (4.40), only the 
coefficients p, Sn a nd S1 2 chan ge. Thi s is be cause wc can determine q, r, s, 531, S32, and S33 independent of p, Sn 
and S 12 from Eqs. ( ^78|) , (JyJ), ( p^) and (|43l|) - (|4~33| ). 



V. THE TOTAL ENERGY AND ANGULAR MOMENTUM 

We take the total energy up to 0(R~ 3 ) at Newtonian order and up to 0(R~ 4 ) at 1PN order. Also the total angular 
momentum up to 0(1) at Newtonian order and up to 0(R~ 1 ) at 1PN order, except for £7dr- 



A. Total energy 



The total energy of the binary system is defined as |2 
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^The change of the physical values is less than 1%, even near the ISCO. 



23 



E 



v 2 U 



DR 



\M 2 R NL 



R 



4 

-0(R- 3 ) 



I J k "I" ^22*) — -Fa-fll* + T^ll* 



(5.7) 



2 

35 



M*ir P nl R FZalal \(-^ - lW 

L v a 2* 7 



1 Ma 



R 2 
' 10 



n 



v*f3 2 



DR 



M*7rpriDR(7A * + a|M2*), 
7M?.R M* 



4 

-0(iT 3 



•^U 2 2 

E V 2 P / P = —M*R f2 DR 
op 



( — 14F a In* + 2F a I 2 2* + I22* 



4^% , 



8M< 



TTjOAo* 



14i? 2 M* 



Egj2^ij2 — —5 



M 3 M 2 



0(R~ 6 ) 



(5.8) 

(5.9) 
(5.10) 
(5.11) 
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In the above equations, we have used the conserved mass 
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instead of the Newtonian mass M, and we have defined R = R/a\. Also we substitute the mean radius of the star 
defined by the conserved mass as 
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and cii* defined by 01*02*03* = a 3 , a2*/ai* = 02/01 and 03*/ai* = 03/01. In Eq. (5.4), 81%, 5Iu, and oW denote 
the perturbed dipole moment, the perturbed quadrupole moment, and the perturbed self-gravity potential energy, 
respectively. These terms are given in Chandrasekhar's textbook [p9j and in our previous paper, (l^] We represent 
these terms in Appendix p[ 



The last term on the right-hand side of Eq. (5.4) is calculated from 
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where the first term on the right-hand side of the above equation denotes the gravitational potential energy which 
the nondisplaced element of star 1 experiences due to the displaced potential of star 2. On the other hand, the 
second term denotes the gravitational potential energy which the displaced element of st ar 1 exp erien ces due to the 
nondisplaced potenti al of star 2. For the Lagrangian displacement vectors given in Eqs. (3.4) - (3.11), the first and 
second terms of Eq. ( 5.16 ) are equal, and calculated as 
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B. Total angular momentum 

The total angular momentum is written as |25| 
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VI. NUMERICAL RESULTS 

In the following subsections, we calculate the equilibrium sequence of the irrotational binary system. 



A. Normalization 



First of all, we introduce non-dimensional parameters according to 
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Using these parameters with the condition So = 0, we can rewrite two groups of equations into non-dimensional forms. 
One of these is ( pjgj ), Q - Q, (|3|) - (g3§), and the other is ( |2~79| ), ( |OoD - ( gig ), fl4~37| ) and ( jQjg ). From 
these fifteen equations, we can determine fifteen variables, p, q, f, s, Sn, S12, S31, S32 and 6*33, and e, /, g, h, S21 
and 5*22. 

After determinat ion of the varia bles, we can calculate the total energy and total angular momentum of the binary 
system from Eqs. (5.1) and (5.18). In the numerical calculations, these values must be normalized. The normalized 
energy and angular momentum are written as 
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where C s is the compactness parameter defined by 

C s 



(6.8) 



We consider the case C s < 1 because of the 1PN approximation. 

At 1PN order, the center of mass of the binary system deviates from the value defined at Newtonian order. The 
center of mass at 1PN order is defined by using the conserved mass as 
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Then, the orbital separation of the 1PN order is written as 



R* = R 



(6.9) 



(6.10) 



(6.11) 



Once an equilibrium sequence is obtained, we search for the minimum point of the total energy. If we find it, we 
call it the ISCO. 



B. Ellipsoidal approximation 



The ellipsoidal approximation, in which the equilibrium configuration is assumed to have an ellipsoidal figure, is 
useful when we study features of stars or binary systems in the 1PN approximation, |l7],[l9|,[54| because the whole 
calculation is carried out by setting Sy = 0, and hence is greatly simplified. The ellipsoidal approximation gives an 
exact solution for a rotating incompressible star in the Newtonian theory. However, it gives only an approximate 
solution for the 1PN case. If the ellipsoidal approximation is truly a robust approximation, it becomes a useful 
method for the study of the 1PN effects. This is the motivation for our investigation of the validity of the ellipsoidal 
approximation. 

In the ellipsoidal approximation, a solution is obtained by setti ng Sg = in a ll th e eq uations. A fter w e set Sg = 0, 
we ca n calculate the velocity potential at 1PN order from Eqs. ( 2.78 ) and (4/7) - (fi~9|), and also ( [2.79 ) and ( 4.10 ) 
(4.12). We note that in the ellipsoidal approximation, the boundary condition (p~5|) is not satisfied. 



C. Results 



The results are shown in Figs. 2-5. In Fig. 2, we represent the total energy and total angular momentum 
as functions of the normalized separation of the binary system i?*/a» and the normalized orbital angular velocity 

\J M*/a%. We can see from this figure that the minimum point of the total angular momentum is slightly different 
from that of the total energy. The deviation between the location of the minimums of the energy and angular 
momentum comes from the fact that we assume a»/i?» is a small parameter and expand the energy up to 0(R~ ) 
and the angular momentum up to 0{R~ r ). In any case, we may expect that the ISCO is located near the two 
minimums. It is found that even if we increase the compactness parameter C s , the angular velocity at the ISCO in 
units of a/ M*/a» has almost the same value as in the Newtonian case, while the orbital separation decreases. This 
feature of the angular velocity at the ISCO is different from that in the case of the corotational binary system. In 
such a case, we have obtained the result that when we increase the compactness parameter, the angular velocity at 
the ISCO increases. ]l9| ] In a previous paper investigating irrotational and incompressible stars, Jig ] we found the 
instability driven by the deformation of 1PN order. We regard such an instability point as the bifurcation point to 
a new sequence. However, in the present study, we do not find such an instability point throughout the equilibrium 
sequence of the binary system until the ISCO. 
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FIG. 2. The total energy and the total angular momentum of the full calculation as functions of the orbital separation 
i?»/a* and the angular velocity fi/(M* /a*) 1 / 2 . The solid, dotted, dashed, and long-dashed lines denote the cases of C s = 
(Newtonian), 0.01, 0.03, and 0.05, respectively. 




FIG. 3. The total energy and the total angular momentum of the ellipsoidal approximation as functions of the orbital 
separation R*/a* and the angular velocity f2/(M* /al) 1 ^ 2 . The identifications of the lines are the same as in Fig. 2. 

In Fig. 3, the total energy and total angular momentum in the ellipsoidal approximation are shown as functions 
of the normalized separation of the binary system and the normalized orbital angular velocity. It is found that the 
quantitative features of the energy and angular momentum are the same as in Fig. 2. Also, we find that values at 1PN 
order (JfipNc, -Bpnc and Jpnc) differ by only a few percent outside the ISCO from those of the full calculation, even 
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if the figure is assumed to be an ellipsoid. This implies that the ellipsoidal approximation gives good approximate 
results. 




FIG. 4. The total energy and the total angular momentum of the approximation in which we neglect the velocity potential 
at 1PN order as functions of the orbital separation i?*/a, and the angular velocity fl/ '(M, / a*) 1 / 2 . The identifications of the 
lines are the same as in Fig. 2. 




FIG. 5. The relation between the angular velocity at the ISCO Qisco(M*/c 3 ) and the compactness parameter C s = M,/c 2 a*. 
The thick and thin lines denote the cases of the full calculation and the absence of the 1PN velocity potential, respectively. 
Filled circles and triangles denote the results of our calculation. Solid, dotted and dashed lines are the fitting lines for n — 
(incompressible case), 0.5 and 0.7, where n denotes the polytropic index. 



29 



However, neglect of the velocity potential of 1PN order has an effect on the quantitative nature of the energy 
and angular momentum. We can see from Fig. 4 that the angular velocity at the minimum point of the energy 
and angular momentum increases when we increase the com pact nes s pa ramet er. F igu re 4 is obtained by se tti ng th e 



veloci ty pot entia l at 1PN order to zero and solving only Eqs. (|4.7|) - (|4.9|) and (|4.3l|) - (|4.35|) , and also ( |4.10|) - ( |4.12| ) 



(4.37) and (4.38) 



In Fig. 5, we show the relation between the compactness parameter C s — Af,/c a, and the angular velocity at the 
ISCO in units of c 3 /M*. Although the 1PN approximation is not valid, we take C s up to 0.25, which we regard as the 
realistic compactness parameter of a neutron star, in order to compare our results with those of Lombardi, Rasio and 
Shapiro |l7j (see Fig. 4 in their paper). In Fig. 5, the thick and thin solid lines denote the cases of the full calculation 
and the absence of the 1PN velocity potential. These lines come from our calculation for an incompressible fluid. The 
dotted and dashed lines denote the "compressible" cases, whose polytropic indices are 0.5 and 0.7, respectively. We 
plot these lines by assuming that the 1PN correction to the angular velocity is the same value as in the incompressible 
case if the stars in the binary system composed of compressible fluid have a polytropic index less than 1, i.e., the hard 
equation of state. Also, we adopt the angular velocity at the ISCO given by Uryu and Eriguchi as that of the 
irrotational and compressible binary system in the Newtonian theory. 

We can see from Fig. 5 that the angular velocity at the ISCO is overestimated when we neglect the 1PN velocity 



potential. Figure 5 agrees qualitatively with that of Lombardi, Rasio and Shapiro. 17 However, our results suggest 
that the angular velocities given by Lombardi, Rasio and Shapiro are overestimated because of the absence of the 
1PN velocity potential in their calculation. The reason that the total force of the system balances with the smaller 
orbital angular velocity than that in the case of the absence of the 1PN velocity potential is as follows. We regard 
the energy minimum point as the ISCO in this paper. As the stars in the binary system become closer and closer to 
each other, the energy terms with a minus sign, such as the binding energy term, decrease and those with a plus sign, 
such as the rotation energy term, increase. At the ISCO, the increase of the energy terms with a plus sign is equal to 
the decrease of those with a minus sign. When the star possesses internal motion, such motion can produce a part of 
the energy with a plus sign. Therefore, if the star with the 1PN velocity potential has the same ellipsoidal figure as 
that without the 1PN velocity potential, the energy of the orbital motion, that is to say, the orbital angular velocity, 
can decrease. 

Finally, we summarize the results in Tables I - VI. In Table I, the angular velocity, energy and angular momentum 
of the Newtonian and 1PN orders are shown along the sequence of the binary system. In this table, the symbol f 
denotes the location of the ISCO of the Newtonian binary system, i.e., the point at which the total energy has its 
minimum. 

In Table II, the total angular velocity, total energy, and total angular momentum at 1PN order are shown. The 
symbol f denotes the location of the ISCO. In Table III, we give the coefficients of the velocity potential and defor- 
mation associated with the star itself, and in Table IV those associated with the binary motion are represented. Also, 
we show the results of the ellipsoidal approximation in Tables V and VI. We can see from Table V that the values of 
the energy, angular momentum and angular velocity at 1PN order in the ellipsoidal approximation deviate by only a 
few percent from those in the full calculation (Table I). 

TABLE I. Equilibrium sequences of the irrotational Darwin-Riemann ellipsoids at Newtonian order and the angular velocity, 
the energy, and the angular momentum at 1PN order along the sequences. The symbol t denotes the energy minimum point 
of the Newtonian sequence, i.e., the ISCO defined in this paper. 



R/ai 


R*/a* 


ai/ai 


az/ai 


SJdr 


En 


Jn 


SO. 


-Bpn 


JpN 


6.00 


6.071 


0.9824 


0.9828 


9.461(-2) 


-1.282 


1.743 


0.1249 


2.944(-2) 


-0.7523 


5.50 


5.584 


0.9771 


0.9777 


0.1073 


-1.289 


1.673 


0.1401 


5.964(-2) 


-0.4832 


5.00 


5.103 


0.9693 


0.9705 


0.1229 


-1.298 


1.600 


0.1584 


9.933(-2) 


-0.1947 


4.50 


4.628 


0.9576 


0.9597 


0.1424 


-1.308 


1.525 


0.1806 


0.1523 


0.1157 


4.00 


4.165 


0.9391 


0.9430 


0.1671 


-1.319 


1.450 


0.2079 


0.2239 


0.4498 


3.50 


3.722 


0.9078 


0.9160 


0.1986 


-1.333 


1.377 


0.2412 


0.3212 


0.8081 


3.00 


3.317 


0.8514 


0.8692 


0.2382 


-1.346 


1.315 


0.2802 


0.4511 


1.190 


t2.582 


3.037 


0.7671 


0.8013 


0.2763 


-1.352 


1.291 


0.3147 


0.5879 


1.548 


2.50 


2.992 


0.7446 


0.7831 


0.2839 


-1.352 


1.292 


0.3215 


0.6225 


1.641 


2.00 


2.842 


0.5586 


0.6238 


0.3226 


-1.319 


1.397 


0.3984 


0.9807 


2.856 
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TABLE II. Equilibrium sequences of the irrotational Darwin-Riemann ellipsoids at 1PN order with the compactness param- 
eters M,jc 2 a, = 0.01, 0.03, and 0.05. The symbol \ denotes the point of the ISCO. 



R/ai 




0,2 J €L\ 


&3 1 0,\ 






E 


J 








M*jc a* = 


U.01 








6.00 


5.987 


0.9824 


0.9828 




9.585(-2J 


-1.282 


1.736 


5.50 


5.506 


0.9771 


0.9777 




0.1087 


-1.289 


1.668 


5 nn 

O.KJKJ 


^ n^n 

O.KJOKJ 


n Qficn 

u.yuyo 


KJ.iJ ( KJO 




0.1244 


1 9Q7 


1 TOR 


4.50 


4.561 


0.9576 


0.9597 




0.1442 


-1.306 


1.526 


4.00 


4.104 


0.9391 


0.9430 




0.1691 


-1.317 


1.454 


3.50 


3.666 


0.9078 


0.9160 




0.2010 


-1.329 


1.385 


3.00 


3.265 


0.8514 


8692 




0.2410 


-1.341 


1.327 


-(-2.618 


3.008 


0.7762 


0.8086 




0.2761 


-1.346 


1.306 


2.50 


2.944 


0.7446 


0.7831 




0.2871 


-1.346 


1.308 


2.00 


2.796 


5586 


6238 




3266 


-1.309 


1.426 








71 /f 1 „ 

M*/c a* = 


n no 
0.03 








6.00 


5.819 


0.9824 


0.9828 




9.835(-2) 


-1.281 


1.721 


5.50 


5.350 


0.9771 


0.9777 




0.1115 


-1.288 


1.658 


5 nn 

O.KJKJ 


*±.ooo 


n Qficn 

u.yuyo 


KJ.iJ ( KJO 




n 1 97R 






4.50 


4.428 


0.9576 


0.9597 




0.1478 


-1.303 


1.529 


4.00 


3.981 


0.9391 


0.9430 




0.1733 


-1.313 


1.463 


3.50 


3.553 


0.9078 


0.9160 




2058 


-1.323 


1.401 


3.00 


3.161 


0.8514 


0.8692 




0.2466 


-1.332 


1.351 


t2.698 


2.958 


0.7951 


0.8238 




0.2747 


-1.335 


1.336 


2.50 


2.848 


0.7446 


0.7831 




2936 


-1.333 


1.341 


2.00 


2.703 


0.5586 


0.6238 




0.3346 


-1.290 


1.483 








71 /f / ^ .-, 

M*/c a* = 


0.05 








6.00 


5.651 


0.9824 


0.9828 




0.1009 


-1.281 


1.706 


5.50 


5.193 


0.9771 


0.9777 




0.1143 


-1.286 


1.648 


5.00 


4.740 


0.9693 


0.9705 




0.1308 


-1.293 


1.590 


4.50 


4.294 


0.9576 


0.9597 




0.1514 


-1.300 


1.531 


4.00 


3.858 


0.9391 


0.9430 




0.1775 


-1.308 


1.472 


3.50 


3.440 


0.9078 


0.9160 




0.2106 


-1.316 


1.417 


3.00 


3.057 


0.8514 


0.8692 




0.2522 


-1.323 


1.374 


t2.795 


2.919 


0.8156 


0.8403 




0.2715 


-1.324 


1.366 


2.50 


2.751 


0.7446 


0.7831 




0.3000 


-1.321 


1.373 


2.00 


2.610 


0.5586 


0.6238 




0.3425 


-1.270 


1.540 



TABLE III. Coefficients of the 1PN velocity potential and the Lagrangian displacement vectors associated with the star 
itself shown along the equilibrium sequence of the irrotational Darwin-Riemann ellipsoid. 



R/ai 


P 


q 


f 




s 




Su 




Sl2 




5*31 




5*32 




S33 


6.00 


-1.634(-2) 


-6.281(-4) 


5.725( 


5) 


4.333( 


5) 


3.205(- 


3) 


-1.997(- 


2) 


-1.950(- 


4) 


-1.851(- 


5) 


6.821(-5) 


5.50 


-1.937(-2) 


-9.253(-4) 


8.352( 


-5) 


5.711( 


-5) 


8.525(- 


4) 


-1.978(- 


2) 


-3.235(- 


4) 


-3.191(- 


5) 


1.152(-4) 


5.00 


-2.296(-2) 


-1.412(-3) 


1.256( 


-4) 


7.251( 


■5) 


-2.921(- 


3) 


-1.884(- 


2) 


-5.599(- 


4) 


-5.826(- 


5) 


2.046(-4) 


4.50 


-2.686(-2) 


-2.248(-3) 


1.948( 


-4) 


8.096( 


•5) 


-8.970(- 


3) 


-1.665(- 


2) 


-1.016(- 


3) 


-1.141(- 


4) 


3.862(-4) 


4.00 


-2.990(-2) 


-3.760(-3) 


3.104( 


4) 


4.680( 


5) 


-1.855(- 


2) 


-1.258(- 


2) 


-1.943(- 


3) 


-2.436(- 


4) 


7.863(-4) 


3.50 


-2.752(-2) 


-6.660(-3) 


4.937( 


4) 


-1.711( 


-4) 


-3.264(- 


2) 


-6.665(- 


3) 


-3.896(- 


3) 


-5.781(- 


4) 


1.760(-3) 


3.00 


-1.683(-3) 


-1.249(-2) 


6.546( 


-4) 


-1.149( 


3) 


-4.573(- 


2) 


-5.496(- 


3) 


-7.864(- 


3) 


-1.548(- 


3) 


4.448(-3) 


2.50 


0.1209 


-2.332(-2) 


-7.605( 


-4) 


-5.280( 


3) 


1 . 1 1 1 (- 


4) 


-5.723(- 


2) 


-1.216(- 


2) 


-4.555(- 


3) 


1.313(-2) 


2.00 


0.6365 


-2.618(-2) 


-2.067( 


-2) 


-1.992( 


2) 


7.509(- 


1) 


-5.522(- 


1) 


2.515(- 


2) 


-1.447(- 


2) 


5.382(-2) 
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TABLE IV. Coefficients of the 1PN velocity potential and the Lagrangian displacement vectors associated with the binary 
motion shown along the equilibrium sequence of the irrotational Darwin-Riemann ellipsoid. 



R/ai 


e 


/ 


9 


h 


S'21 


S22 


6.00 


0.1750 


-2.520(-2) 


-3.005(-2) 


-2.997(-2) 


1.465(-3) 


-1.608(-3) 


5.50 


0.2113 


-2.517(-2) 


-3.178(-2) 


-3.164(-2) 


2.264(-3) 


-2.484(-3) 


5.00 


0.2567 


-2.465(-2) 


-3.394(-2) 


-3.368(-2) 


3.649(-3) 


-3.999(-3) 


4.50 


0.3142 


-2.320(-2) 


-3.678(-2) 


-3.627(-2) 


6.185(-3) 


-6.766(-3) 


4.00 


0.3880 


-1.991(-2) 


-4.077(-2) 


-3.970(-2) 


1.116(-2) 


-1.216(-2) 


3.50 


0.4835 


-1.278(-2) 


-4.697(-2) 


-4.450(-2) 


2.175(-2) 


-2.351(-2) 


3.00 


0.6057 


3.204(-3) 


-5.798(-2) 


-5.167(-2) 


4.663(-2) 


-4.944(-2) 


2.50 


0.7498 


4.035(-2) 


-8.037(-2) 


-6.307(-2) 


0.1099 


-0.1105 


2.00 


0.8687 


0.1120 


-0.1227 


-8.151(-2) 


0.2634 


-0.2280 



TABLE V. The 1PN angular velocity, energy, and angular momentum in the ellipsoidal approximation along the equilibrium 
sequence of the irrotational Darwin-Riemann ellipsoid. The subscript "e" denotes the case of the ellipsoidal approximation. 



R/ai 




R*/a, 






sn c 




-BpNe 




JpNc 


6.00 




6.071 






0.1249 




2.940(-2) 




-0.7527 


5.50 




5.584 






0.1401 




5.961(-2) 




-0.4835 


5.00 




5.103 






0.1584 




9.935(-2) 




-0.1945 


4.50 




4.628 






0.1808 




0.1524 




0.1167 


4.00 




4.165 






0.2084 




0.2244 




0.4526 


3.50 




3.722 






0.2426 




0.3225 




0.8149 


3.00 




3.317 






0.2837 




0.4541 




1.204 


2.50 




2.992 






0.3276 




0.6165 




1.626 


2.00 




2.842 






0.3527 




0.7600 




2.153 


TABLE VI. Coefficients of the 1PN velocity potential 


in the ellipsoidal approximat 


ion along the equilibrium sequence of the 


irrotational Darwin-Riemann ellipsoid. 


The subscript "e 


' denotes the case of the ellipsoidal approximation. 




R/ai 


Pc 


?o 




fc 


So 


e c 


fc 


9c 




6.00 


-1.853(-2) 


-6.155(-4) 


4.426(-5) 


4.445(-5) 


0.1750 


-2.537(-2) 


-3.000(-2) 


-2.997(-2) 


5.50 


-2.157(-2) 


-9.015(-4) 


5.892(-5) 


5.939(-5) 


0.2113 


-2.546(-2) 


-3.168(-2) 


-3.164(-2) 


5.00 


-2.489(-2) 


-1.365(-3) 


7.629(-5) 


7.751(-5) 


0.2567 


-2.519(-2) 


-3.376(-2) 


-3.368(-2) 


4.50 


-2.791(-2) 


-2.147(-3) 


8.957(-5) 


9.302(-5) 


0.3141 


-2.427(-2) 


-3.642(-2) 


-3.628(-2) 


4.00 


-2.881(-2) 


-3.529(-3) 


6.887(-5) 


7.960(-5) 


0.3879 


-2.220(-2) 


-4.001(-2) 


-3.972(-2) 


3.50 


-2.218(-2) 


-6.096(-3) 


-1.043(-4) 


-6.787(-5) 


0.4833 


-1.810(-2) 


-4.518(-2) 


-4.455(-2) 


3.00 


8.177(-3) 


-1.107(-2) 


-8.948(-4) 


-7.645(-4) 


0.6055 


-1.052(-2) 


-5.335(-2) 


-5.185(-2) 


2.50 


0.1068 


-2.057(-2) 


-4.058(-3) 


-3.635(-3) 


0.7515 


2.764(-3) 


-6.759(-2) 


-6.380(-2) 


2.00 


0.3336 


-3.347(-2) 


-1.531(-2) 


-1.414(-2) 


0.8843 


2.684(-2) 


-9.351(-2) 


-8.377(-2) 
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VII. SUMMARY AND DISCUSSION 



A. Summary 

Using the scheme developed in a previous paper, |Q we have investigated the equilibrium sequences of irrotational 
and incompressible binary systems in the 1PN approximation. Our results presented in the previous section should 
be useful to check the numerical code for solving the irrotational binary problem. The conclusions are as follows. 

(1) Due to the 1PN effect, the orbital separation at the ISCO decreases in proportion to the compactness parameter 
Af*/c 2 a». 

(2) The orbital angular velocity at the ISCO in units of yj M*/a| has almost the same value as in the Newtonian case, 
even if we increase the compactness parameter. 

(3) The ellipsoidal approximation gives fairly accurate results throughout the equilibrium sequence. 

(4) It is important to include the velocity potential at 1PN order. 



B. Discussion 



Finally, we discuss the velocity potential at 1PN order. As stated in a previous paper, |Q the internal velocity 
has an X3 component. Therefore, we suggest that when we solve an irrotational binary problem numerically, we must 
include the X3 component of the internal velocity. 

Next, we discuss the effect of the deformation at Newtonian order. Even in the incompressible case, stars in the 
binary system deform from ellipsoidal figures at Newtonian order. |l3| We can estimate the contribution from the 
octupole moment on the angular velocity. First, we assume that the deformation from an ellipsoidal figure is a small 
perturbation. Accordingly, the deformation is expressed by the Lagrangian displacement vectors (3.10) and (3.11) 
for the octupole moment. Here we neglect the contribution of the hexadecapole and other moments as higher order 
terms. Then, the angular velocity is written 



9. 2 



I8J11 
R 5 



MR 



2M 

4 

~R? 

+\sn 2 + o(r- j ), 



2M 



R 5 



(2 — F a )_F a fi DR 



(SI 



1 Jn 



[2(Miii) n - 3(« 122 ) N - 3((SIi 33 ) 



(7.1) 



where the subscript N denotes the perturbation of Newtonian order. Equation (7.1) is approximated as 



n 2 = 



2M 



5R 2 



2a 2 a 3 



(7.2) 



where ti\ and i 22 denote the amplitude of the Lagrangian displacement vectors ( 3.10 ) and ( 3.11 ), and 



w 2 



1 

5R 


1 -f 


3 

" 5R 2 


_l_ 


6 


('- 


35.R 3 


12 


(al- 


35R 3 



<4) 



3 «> \ 



(7.3) 
(7.4) 



From Table II, we can obtain the values in the range of 0.01 < C s < 0.05 as 



w 2 i 22 



0.0243 < a 



0.14*21, 
0.001? 22 , 
R 3 



' 2a 2 a3 



■sn 2 < 0.121. 



(7.5) 
(7.6) 

(7.7) 



Therefore, if the amplitude ti\ (and i 22 ) is much less than 1, i.e., in the range of < i 2 i < 0.1, the contribution of 
the octupole deformation to the angular velocity is smaller than that of the 1PN correction. On the other hand, if t 2 i 
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and/or £22 have values near or greater than 1, the above perturbative treatments break down. In this case, we cannot 
estimate the contribution of the octupole moment until we construct the equilibrium figure including the octupole 
moment. 

Moreover, we would like to mention the importance of the velocity potential of 1PN order. If we neglect it, we 
obtain different features for the ISCO. Accordingly, we must take into account the velocity potential at 1PN order 
even when we use the ellipsoidal approximation. 

From Fig. 5, we can estimate the frequency of gravitational waves at the ISCO as 



450 [Hzl <, 



1.6M 



few £ 800 [Hz], 



(7.; 



if we consider that we have 55 n 5s 1 for the polytropic indicies and 0.15 Ss C s Ss 0.25 for the compactness parameter 
of neutron stars. 
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APPENDIX A: EXPLICIT FORM OF SU 



The correction of the self-gravity potential induced by the deformation of the binary system is given by 



(Al) 



where 



suw = -yJL f d 



\x — x'\ 



(A2) 



Substituting the Lagrangian displacement vectors into Eq. (A2), we separately obtain SU as 

SU = SU 1 ^ 1 + 5U 2 ^\ 



(A3) 



where 



SU 1 



SU 



>S'ii(-D3,3 — + Si2(D 3 3 — 1)2,2) + Ssx ^112,2 — 2^111,1 

+S32 (-D223.3 — -£*222,2^ + S33 ^-Dl33,l — -1^333,3 

"Sotffr 1 + ^21(^13,3 - l D n,i) + S 22 {D 1 ^ - D 12>2 ) 
"SiiG ^ 1 - d ia V ) + Si2{Dl^ 1 - D^ 1 ) + S3i(Diu% - ^iim 

+ •532 ^223,3 _ 2-^222,2) + ^33^-0133^1 _ 2^333^3 



(A4) 



(A5) 



Note that the sign of & for <5£/ 2_+1 is opposite to that of SU 1 ^ 1 if the components of & are even functions of Xj 
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APPENDIX B: PERTURBED TERMS 

The definitions of Sli, Slij, and SW are 



where 



The components are 



6Wu 



SW- 



22 



SW- 



33 



SI, 

Shi = 



/ 



i) 



SW = J2SW U , 

i=l 

JWij = np(-2B l: jSI l: j + a?6ij A U SI U ). 



5h = j^(MS + S 2 iIu), 
Sin = ^2 ^n/n + -531/iin — S'33/1133^ , 

<^22 = ^2 (^12^22 — 5 , 3i/n22 + g 532/2222^) J 

^2 ^ — 5n/33 — ^12/33 — 532/2233 + ^533/3333 



C 
27T/9 



27T/9 



27T/9 



5n{(a^u - 2Bii)/n - ^13/33} 

+5 i2 (af A12/22 - af^is^ss) 

+5 3 i{ ^(af^n - 2 J B n )/ im - 0^12/1122} 

+532 (g a l ^12/2222 ^ 0^13/2233) 

+5 33 {-(a^ii - 2Bn)/ii33 + ia?Ai3/ 3 33 3 | 

5n(a2Ai 2 /ll - ft2^23^33) 

+ 5l2|(a2A 2 2 — 2/?22)/22 — 0^23/33 1 

+ 531 j -02^12^1111 — {a\A 2 2 — 2-B 2 2)/ll22| 

+532| g(a 2 ^22 — 2_B 2 2)/2222 — a!^23-^2233 j 

+533^-02^12/1133 + ga2^23^3333) j 

5u {0^13/11 - (0^33 - 2B33)/33} 
+512(03^23/22 - (a^ 33 - 2 J B 3 3)/3 3 } 
+531^-03^13/1111 — 03^23/1122) 

+5 32 { — 03^23/2222 _ (O3A33 — 2/?33)/2233 
+5 33 { -03^13/1133 + -(03^33 - 2Z?3 3 )/3333 j 
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